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Abstract 

Suppose F and G are two life distribution functions. It is said that F is more IFRA than 
G (written by F <* G) if G~ 1 F(x) is starshaped on (0, oo). In this paper, the problem of 
testing H 0 : F =* G against Hi : F <* G and F 7 ^* G is considered in both cases when 
G is known and when G is unknown. We propose a new test based on U-statistics and 
obtain the asymptotic distribution of the test statistics. The new test is compared with 
some well known tests in the literature. In addition, we apply our test to a real data set 
in the context of reliability. 

Keywords : Asymptotic normality, star order, increasing failure rate average, Pitman’s 


asymptotic efficiency, U-statistic. 



1 Introduction 


Let A" be a lifetime of an appliance with density function /, distribution function F and 
survival function F. Let also denote F~ x as the right continuous inverse function of F. X 
is said to be IFRA (increasing failure rate average) if f F (x) = — l f in '" is nondecreasing in 
x > 0 which is equivalent to that — log ^b) j g nondecreasing in x > 0 where r F (x) = It 
is of considerable interest to producers and users of the appliances to evaluate the severity 
of average failure risk at a particular point of time and to see if r F (x ) is either increasing 
or decreasing in time. That is, it is of practical importance to characterize the aging class 
of underlying random lifetimes. In particular, since the IFRA class of aging is one of 
the most important aging classes, testing that the distribution F has a constant hazard 
rate against the hypothesis that F is IFRA has been studied extensively in the literature; 
see for example, Deshpande (1983), Kochar (1985), Link (1989), Ahmad (2000) and El- 
Bassiouny (2003) among others. In fact, F is IFRA if and only if Ex is nondecreasing 
in x > 0 or equivalently is nondecreasing in u G (0,1) where E\ is an exponential 

distribution with mean A. This implies that F ages faster than E , i.e., F is more IFRA 
than E\. 

In order to evaluate the performance of an appliance, we need to compare its aging 
behavior with some distributions other than exponential distribution such as the Weibull, 
gamma, linear failure rate or even an unknown distribution G. The notion of the star 
order that establishes an equivalent class of distributions is one of the useful tools for this 
comparison. Let Y be another non-negative random variable with distribution function 
G. We say that X is less than Y with respect to the star order (written by X <* Y or 
F <* G ) if G~ 1 F(x ) is starshaped on [0,oo); that is, G * — is nondecreasing in x > 0. 
It is known that 

MF- 1 ^)) ; 


F <* G -v=> 


is nondecreasing in u G (0,1), 


( 1 . 1 ) 


fciG-^u)) 

where r F and fa are failure rate average functions of F and G, respectively. Using (II.IF 
the relation A" <* Y is interpreted as X ages faster than Y and it is said that A" is more 
IFRA than Y (cf. Kochar and Xu, 2011 ). It is obvious that if F <* G and G <* F then 
F(x) = G(ax) for all x > 0 and some a > 0. In this case, we say F =* G. 
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Izadi and Khaledi (2012) have considered the problem of testing the null hypothesis 
H 0 : F =* G against II\ : F <* G and F 7 A G. They proposed a test based on kernel 
density estimation. In this paper, we further study this problem of testing in the one- 
sample as well as the two-sample problem and propose a new simple test based on a 
U-statistic. In both cases, we compare the new proposed test with some well known tests 
in the literature. It is found that our test is comparable to the others. 

To establish our new test we need the following lemma. 


Lemma 1.1 Let Xi,X 2 (Y \, Y 2 ) be two independent copies of the random variable X (Y) 
with distribution function F (G) and let /ip = E[max{Xi, X 2 }] (p G = ^[maxjYi, Y 2 }]) 
where E[.] is the expectation operator. If F is more IFRA than G, then 


( 2 ) ( 2 ) 
Of < Eg 

Of Og 

where p F (og) is the expectation of F (G). 


Proof: We know that more IFRA order is scale invariant. Thus, X <* Y implies <* 
Y. Now, the required result follows from Theorem 7.6 of Barlow and Proschan (1981, page 
122 ). 

Remark 1.1 The above lemma has been proved by Xie and Lai (1996) under the condition 
that F is more IFR than G (for definition, see Shaked and Shantikumar, 2007, p. 21 f) 
which is stronger than more IFRA order. 

,P) // (2) 

Now, let S F = S G = and 

5(F,G) = 5 f - 5 g . (1.2) 

It is obvious that if F =* G, then 8(F, G) = 0 and if F <* G and F 7 ^* G, then it follows 
from Lemma 11.11 that S(F,G) < 0. That is, S(F,G) can be considered as a measure of 
departure from H 0 : F =* G in favor of Hi : F <* G and F 7 ^* G. So, our test statistic is 
based on the estimation of S(F, G). 

The organization of this paper is as follows. In Section 2, we propose the new test for 
the case when G is known. The case when G is unknown is studied in Section 3. In Section 
4, the performance of our test is evaluated and compared. 
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2 The One-Sample Problem 

Let Go be a known distribution function and Xi,..., X n be a random sample from an 
unknown distribution F. Now by using the measure (11.21) . the test statistic 

&Go 


is used for testing 


H 0 :F =* Go 


against 


Hi : F Gq and F 7 G Go 


where 


$f — 


EE max {Ad, Xj} 

*7G 


(2,3) 


n(n — 1)X 

and X is the mean of the random sample. In the next theorem, we obtain the asymptotic 
distribution of S(F, G 0 ) by using the standard theory of U-statistics. 

Theorem 2.1 Suppose E[max{X 1: X 2 } — ^-(X 1 + X 2 )} 2 < 00. As n —> 00, n ly/2 [<5(F, G 0 ) — 
S(F, G 0 )] is asymptotically normal with mean 0 and variance 

1 / A \ 

a\ = — xVar ( XF(X) + / tdF(t) - -f A j . (2.4) 

l l F \ Jx 2 J 

Under H 0 , n 1//2 5(F, Go) is asymptotically normal with mean 0 and variance ctq = Oq . 


Proof: First note that 


5f — 5f — 


EE 


max{A i ,A i }-y(A i + A j ) 


n(n — 1)X 


iAj 


n(n — 1)X 
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where 


Xj) = max {A",, Xj} - y (X t + X)). 

Let define 

rj~t* _ ^7^3 

n{n — 1) 

By the standard theory of U-statistics, if E[(p 2 (X 1 , X 2 )\ < oo, as n — » oo 

^4 ^- 4 iV(0,1) 

CT * 

where 

4 = 4 x Var(MX)) 

and 

M x ) = ^[ 04 ,^)]- 

Now by the strong law of large numbers we have X 4' /i F and hence, by Slutsky theorem 

/s 2 

yfn[5(F, G 0 ) — 8(F,G 0 )] is asymptotically normal with mean 0 and variance cr F = jp. 
Under H 0 , S(F, G 0 ) = 0 and a g = Gq q . □ 


A small value of 5(F,G o) indicates that testing H 0 against H i is significant. Thus, 
we reject H 0 at level a if n 1//2 <5(F, G 0 )/ctgo < z ai where z a is a th quantile of the standard 
normal distribution. 

In the case Gq(x ) = E\(x) = 1 — exp(— Xx), x > 0 and A > 0, the problem is testing 
the null hypothesis H 0 : F is an exponential distribution against the alternative hypothesis 
H\ : F is IFRA and not exponential. It can be shown that 5 Fx = f and a F =4- By 
the above theorem, under H 0 , ^Jn(5 F — 3/2) is asymptotically normal with mean 0 and 
variance 4- Thus we reject H 0 in favor of H± if \/l2 n(8 F — 3/2) < z a . 

In the following we find the exact distribution of 5 F under the hypothesis F is an 
exponential distribution. First, note that we can rewrite S F as 



2 - W-l 


i=l 

n{n — 1)X 


n 

i= 1 
n 

i —1 


(2.5) 
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where is the i th order statistic of AYs, 


D i=(n-%+ l)(X (i) - c i:n = (n _j)( n _,- + 1 ) 

and assuming A( 0 ) = 0. Now, by the same arguments as in Langenberg and Srinivasan 
(1979), we will get the following result. 


Theorem 2.2 Let F be an exponential distribution, then 


P{S F < x} = 1 - Y\ C ' n 1 I(x < c i:n ) (2.6) 

, , Q:n Cj:n 

1=1 j=l J 

where /(.) is the usual indicator function. 


By using Theorem 12.21 we tabulate the critical point of ^12 n(5p — 3/2) under expo- 
nentiality for small sample sizes (< 40) in Table [B So, for small sample sizes, we reject 
exponentiality in favor of IFRA-ness if \/12 n(5 F — 3/2) is smaller than the critical point 
in Table |T] corresponding with the level of significance chosen. 

El-Bassiouny (2003) has considered the problem of testing exponentiality against IFRA- 
ness in the alternative and proposed a class of test. His test is based on the test statistics 


A 


i<j 


2 V V f min{A7 +1 


A7 +1 } - 


X, 


r+1 


r+1 


n(n - l)A r+1 


and large values of A r+1 are significant for the considered problem of testing. If r = 0, 

2 - y) 

Ai _ »<j - - 

n(n — 1)X 

i<j 1 

n{n — 1)X 2 


5 


(2.7) 











Table 1: 


Critical values of y/12n(5p — 


3/2) for small sample sizes 


n 

0.01 

a: Lower Tail 

0.05 

0.1 

0.1 

a: Upper Tail 
0.05 

0.01 

2 -2.400500 

-2.204541 

-1.959592 

1.959592 

2.204541 

2.400500 

3 -2.575752 

-2.051328 

-1.658360 

1.658360 

2.051328 

2.575752 

4 -2.560006 

-1.918143 

-1.516280 

1.516280 

1.918143 

2.560006 

5 -2.517587 

-1.846175 

-1.458997 

1.458997 

1.846175 

2.517587 

6 -2.482569 

-1.807959 

-1.424446 

1.424446 

1.807959 

2.482569 

7 -2.458901 

-1.781575 

-1.400767 

1.400767 

1.781575 

2.458901 

8 -2.441786 

-1.762473 

-1.383974 

1.383974 

1.762473 

2.441786 

9 -2.428500 

-1.748106 

-1.371312 

1.371312 

1.748106 

2.428500 

10 -2.417939 

-1.736865 

-1.361442 

1.361442 

1.736865 

2.417939 

11 -2.409356 

-1.727862 

-1.353531 

1.353531 

1.727862 

2.409356 

12 -2.402239 

-1.720450 

-1.347047 

1.347047 

1.720450 

2.402239 

13 -2.396243 

-1.714193 

-1.341635 

1.341635 

1.714193 

2.396243 

14 -2.391124 

-1.708937 

-1.337050 

1.337050 

1.708937 

2.391124 

15 -2.386703 

-1.704422 

-1.333116 

1.333116 

1.704422 

2.386703 

16 -2.382846 

-1.700502 

-1.329703 

1.329703 

1.700502 

2.382846 

17 -2.379451 

-1.697066 

-1.326714 

1.326714 

1.697066 

2.379451 

18 -2.376441 

-1.694029 

-1.324074 

1.324074 

1.694029 

2.376441 

19 -2.373754 

-1.691327 

-1.321727 

1.321727 

1.691327 

2.373754 

20 -2.371340 

-1.688906 

-1.319625 

1.319625 

1.688906 

2.371340 

21 -2.369160 

-1.686725 

-1.317732 

1.317732 

1.686725 

2.369160 

22 -2.367182 

-1.684749 

-1.316018 

1.316018 

1.684749 

2.367182 

23 -2.365378 

-1.682952 

-1.314460 

1.314460 

1.682952 

2.365378 

24 -2.363726 

-1.681309 

-1.313036 

1.313036 

1.681309 

2.363726 

25 -2.362209 

-1.679803 

-1.311730 

1.311730 

1.679803 

2.362209 

26 -2.360810 

-1.678415 

-1.310529 

1.310529 

1.678415 

2.360810 

27 -2.359516 

-1.677134 

-1.309419 

1.309419 

1.677134 

2.359516 

28 -2.358316 

-1.675947 

-1.308392 

1.308392 

1.675947 

2.358316 

29 -2.357199 

-1.674844 

-1.307437 

1.307437 

1.674844 

2.357199 

30 -2.356158 

-1.673817 

-1.306548 

1.306548 

1.673817 

2.356158 

31 -2.355185 

-1.672857 

-1.305718 

1.305718 

1.672857 

2.355185 

32 -2.354273 

-1.671960 

-1.304942 

1.304942 

1.671960 

2.354273 

33 -2.353418 

-1.671118 

-1.304214 

1.304214 

1.671118 

2.353417 

34 -2.352613 

-1.670326 

-1.303529 

1.303529 

1.670326 

2.352612 

35 -2.351855 

-1.669581 

-1.302885 

1.302885 

1.669581 

2.351854 

36 -2.351140 

-1.668878 

-1.302278 

1.302278 

1.668878 

2.351138 

37 -2.350454 

-1.668214 

-1.301704 

1.301704 

1.668214 

2.350461 

38 -2.349821 

-1.667586 

-1.301161 

1.301161 

1.667586 

2.349821 

39 -2.349184 

-1.666990 

-1.300647 

1.300647 

1.666991 

2.349213 

40 -2.348688 

-1.666426 

-1.300159 

1.300159 

1.666426 

2.348636 
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On the other hand, using the fact that Sp x = §, 


S(F,E X ) 


-6 


Ex 


Y Y maX { X " X A 
_ 

n(n — 1)X 

[{ X i + X j) — min{Xj, Xj}\ 

j<j _ 

n(n — 1)X 


3 

2 


V ]T(A, + Xj) - ^ ^ 


- i<j 


i<j 


n 


n(n — 1)X 

n 


2—1 


i<j 


n(n — 1)A" 


3 

2 


2^^ min {X i ,X i } 

1 i<j 

2 n(n - 1 )A 

-A,. 


3 

2 


That is, for the case when r = 0 and Go is an exponential distribution, the proposed test 
is equivalent to that of El-bassiouny (2003). 

It is worth to mention that our test is consistent; that is, if /3 n {F) is the power of our 
test, then under the alternative hypothesis, lini^oo / 3 n (F) = 1 which follows from Theorem 
12.11 and Problem 2.3.16 in Lehmann (1999), 


3 The Two-Sample Problem 

In this section, we consider the two-sample problem when G is unknown. Let X \, X 2l ..., X n 
and Yi,Y 2 , , Y m be two independent random samples from unknown distribution func¬ 
tions F and G, respectively, and N = n + m. Assume that 5p is as in (j2.3[) and Sq is 
defined similarly in terms of Yi,..., Y m and Y. The test statistic 

5(F,G) = S f ~5 g 

which is the estimate of the measure in (11.21) is used for testing the null hypothesis 

H 0 :F=*G (3.8) 
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against the alternative hypothesis 


Hi : F <* G and F G. 


(3.9) 


Small values of S(F, G) are significant for testing Ho against H\. In the following theorem 
we obtain the asymptotic distribution of S(F,G). 


Theorem 3.1 If E[ max{Xi, X 2 } - f (X x + X 2 )} 2 and £[inax{Ti, Y 2 } - f (Y 1 + 1" 2 )] 2 are 
finite mid n and m —> oo such that ^ —> c, c G (0, |], i/ien y/N(d(F,G) — S(F,G)) is 
asymptotically normal with mean 0 and variance 


N 


N 


a Fr = —(j F H- a r 

F ’ G n F m G 


where a 2 F is given in \2.f\) and a G is defined similarly in terms of Y 

Proof: It is easy to see that 

y/N(S(F, G) - S(F, G)) y/mla F 


cf,g 


\Jmo 2 F + na G 
\[no G 


{s/nidp — dp)/crp] 
[Vm(S G - d G )/a G \- 


y/ map + na\ 

From the result of Theorem 2.1, as both n and m —> oo, we have that 


\/n(fiF — 5 f)/ctf iV(0,1) and y/m{d G — 5 G )/a G -4- iV(0,1). 

Since dp and d G are independent, the required result follows from the fact that convergence 
in distribution is closed under the convolution of independent sequences of random variables 
(cf. Theorem 6.6 of Gut (2009), page 169). ■ 

In practice crf G is unknown, but it can be estimated by the consistent estimator 


i.2 _N 2 N 2 
ap ’ G ~ ^° F + m a ° 


(3.10) 


where 


Y.Yx.) 


&F — 4 x 


i =1 


nX 2 

n 


fcXi) = + 












and 


E+ y <) 


er, 


G 


= 4 x 


i=1 


mY 2 


m S 

w) = r-^E[ max { y <> r j}-f w+n)i- 


m 


i=i 


Now by Slutsky theorem, under H 0 , \/N5(F,G)/Jf,g is asymptotically normal with mean 
0 and variance 1 as both n and m —> oo. Hence, for large sample sizes, H 0 is rejected at 
level a if ^/N5{F,G)/Jf,g < +*• 


4 Simulation Study 

In this section, we study the performance of our test and compare it with some well known 
tests in the literature for the one-sample and the two-sample problems. 


4.1 The One-Sample 

We recall that in the one-sample problem we consider testing H 0 : F =* Go against 
Hi : F <* Go and F 7 G Go when Go is a known distribution. For the case when Gq(x) = 
1 — exp{—Ax}, x > 0 , we compare our proposed test with the following well known tests 
which are in the literature. Note that in this case the problem is testing exponentiality 
against IFRA-ness. 


Deshpande (1983): The test statistics is 


J b 


1 

n(n — 1 ) 


EE hb (Xj , Xj ), 


be ( 0 , 1 ) 


where hb(x,y ) = 1, if x > by, 0, otherwise. Large values of Jb are used to reject 

exponentiality in favor of IFRA-ness. It has been shown that under H 0 , n l ^ 2 {Jb — 

(b + ly 1 ) is asymptotically normal with mean zero and variance 4£ where 

r _ 1 b 1 2(1 — b) 2b 4 1 

^ “ 4 {1 + V+2 + 26+1 + 6 + 1 “ 6 2 + 6 + 1 ~ (6 + l) 2 ^' 

Deshpande (1983) has recommended 6 = 0.9. 
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Kochar (1985): H 0 is rejected for large values of 

E”=, •/(*)*« 


T = 

-L 71, 


nX 


J(u) = 2(1 — u)[l — log(l — «)] — 1. 


(4.11) 


The asymptotic distribution of (108u/17) 1//2 T n is the standard normal distribution. 


Link (1989): Large values of the test statistic 

r = 2 


EE 


A', 


(0 

n(n - 1) X {j) ' 


(4.12) 


certify that F is 1FR.A. For large values of n, under H 0 , the distribution of ^ 

is approximately standard normal. 


Ahmad (2000): The test statistic is 






(4.13) 


where A; is a known symmetric density function and a n is a sequence of positive real 
numbers such that na n —» oo and naf z —> 0. Under some conditions, — |) 

is asymptotically normal with mean zero and variance 1, when F is an exponential 
distribution (Ahmad, 2000). F[ 0 is rejected at level a if — 1) > z\- a . 

Ahmad has recommended standard normal density as kernel function and a n = n~^. 


First, we investigate the accuracy of normal distribution as the limit distribution of the 
test statistics under H 0 . In order to do this, we simulate the size of the tests for nominal 
sizes a— 0.01, 0.05, 0.1 and large sample sizes n = 40(5)60(10)70. In the simulation, 
10000 samples are generated from exponential distribution with mean 1. The calculated 
size is the proportion of 10000 generated samples that resulted in rejection of H 0 where 
the rejection regions have been obtained by the asymptotical distribution of test statistics. 
The simulated values were tabulated in Table [2j All simulations were done by R package. 

From Table [2J we find that the tests by Deshpande (1983) and Kochar (1985) are 
over shoot the nominal sizes for all sample sizes. The simulated sizes of the tests due to 
Link (1989) and Ahmad (2000) are greater than the nominal sizes but Link’s test always 
dominates Ahmad’s test. It is clear from the contents of Table [2] that the simulated sizes 
of our new test are much closer to the nominal sizes for all sample sizes. 
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Table 2: Simulated sizes of our test for different nominal sizes and large sample sizes 


n 


0.01 

nominal size (a) 
0.05 

0.1 

n 


0.01 

nominal size (a) 
0.05 

0.1 

40 

S(F,E) 

0.0104 

0.0518 

0.1044 

55 

S(F,E) 

0.0108 

0.0500 

0.1003 


Jq.9 

0.0637 

0.1243 

0.1709 


Jo.9 

0.0550 

0.1101 

0.1674 


T n 

0.0396 

0.1815 

0.3157 


T n 

0.0346 

0.1565 

0.2753 


r 

0.0181 

0.0612 

0.1110 


r 

0.0151 

0.0569 

0.1063 


a f 

0.0311 

0.0772 

0.1196 


a f 

0.0324 

0.0783 

0.1237 

45 

S(F,E) 

0.0106 

0.0505 

0.1015 

60 

5(F,E) 

0.0101 

0.0517 

0.1042 


Jo.9 

0.0661 

0.1163 

0.1769 


Jo.9 

0.0494 

0.1128 

0.1680 


T n 

0.0380 

0.1704 

0.2938 


T n 

0.0349 

0.1548 

0.2666 


r 

0.0167 

0.0584 

0.1089 


r 

0.0157 

0.0581 

0.1095 


a f 

0.0338 

0.0796 

0.1232 


a f 

0.0304 

0.0742 

0.1209 

50 

S(F,E ) 

0.0112 

0.0494 

0.1009 

70 

5(F,E) 

0.0090 

0.0489 

0.1024 


Jo.9 

0.0601 

0.1220 

0.1736 


Jo.9 

0.0469 

0.1048 

0.1587 


T 

J-n 

0.0371 

0.1654 

0.2852 


T n 

0.0303 

0.1410 

0.2534 


r 

0.0166 

0.0563 

0.1055 


r 

0.0147 

0.0580 

0.1072 


a f 

0.0312 

0.0810 

0.1232 


a f 

0.0310 

0.0783 

0.1240 


In the following, to assess how our proposed test performs relatively, we first con¬ 
sider the large sample sizes and use the measure of Pitman’s asymptotic relative efficiency 
(PARE) (cf. Nikitin, 1995, Section 1.4). Consider testing H 0 that F is an exponential 
distribution against Hi that F = Fg n where 0 n = 0 o + kn ~^, k is an arbitrary positive 
constant and Fg 0 is exponential. Then, Pitman’s asymptotic efficiency (PAE) of a test 
based on statistic T n is 


PAE(T n ) 


lim 

n—too 


9Eg (T n ) | 

09 10=00 

Va,r eo [-JnT„] 


(4.14) 


Using (14.1411 . the PAE of our test is given by 


PAE(5(F, E\)) = 

We consider three families of Weibull, Linear failure rate and Makeham distributions with 
the following density functions. 


/ 9S f , 

\~dF 10=e 


cr 


Ee 0 


(1) Weibull Distribution: 


fo(x) = 9x e 1 e x °, x > 0, 9 > 1. 
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(2) Linear Failure Rate Distribution: 


Q x ^ 

fg(x) = (1 + 9x)e~ x 2 “, x > 0, 9 > 0. 


(3) Makeham Distribution: 

fe{x) = (1 + 9(1 - e - x )) e -x-<X.x+ e - x - 1 ) > x > o, 9 > 0 . 

PAE of our test ( 8(F,E \)), Deshpande’s test («/&), Kochar’s test (T n ), Link’s test (r) and 
Ahmad’s test (A F ) are presented in Tabled . In Tabled! PARE of our test with respect 
to the others has been obtained. It is observed that our test dominated the others except 
Kochar’s test for the LFR alternative case. 


Table 3: PAE of 5(F,E \), J 0 . 9 , T n , F, A F . 


Test \ Hi 

Weibull 

LFR 

Makeham 

S(F,E X ) 

1.4414 

0.75 

0.0833 

Jo. 9 

1.35 

0.3369 

0.0666 

T 

J- n 

1.247 

0.8933 

0.0784 

r 

1.3867 

0.2681 

0.0563 

A f 

1.35 

0.3375 

0.0667 


Table 4: PARE(h(F, E), T) = T = J 0 . 9 , T n , T, A F . 


Test\ Hi Weibull LFR Makeham 


J0.9 

1.0677 

2.222 

1.2489 

T 

1.1558 

0.8396 

1.0625 

r 

1.0394 

2.7974 

1.479 

a f 

1.0677 

2.222 

1.2489 


In practice, the available samples are small. So, it is important to investigate the 
power of the tests and compare them for small sample sizes. Proportion of 10000 samples 
(with small sizes 5(3)15) that reject exponentiality in favor of IFRA-ness is considered 
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for estimating the power of the tests, fn the alternative, we consider Weibull, LFR and 
Makeham distributions. The critical points of J 09 , T n , T and at signihcance level 
a = 0.05 for small sample sizes have been derived from their corresponding papers. Table 
[5] shows the simulated powers of the tests for different alternatives. It is observed that 
in Weibull and Makeham alternatives our new test is more powerful than the others in 
all sample sizes. In the LFR alternative, Kochar’s test dominates the other tests while 
our proposed test is comparable. Also, Kochar’s test and Link’s test are comparable in 
Weibull and Makeham alternatives and are more powerful than the tests of Deshpande 
and Ahmad. 


Table 5: Simulated power of the tests at level of signihcance 0.05 for small sample sizes. 


n 



1.2 

Weibull(0) 

2 

3 

0.2 

LFR(0) 

1 

2.5 

0.2 

Makeham(0) 

1 

2.5 

5 

S(F, 

Ex) 

0.0902 

0.3886 

0.7496 

0.0645 

0.1029 

0.1422 

0.0579 

0.0821 

0.1111 


Jo. 9 


0.0658 

0.1782 

0.3157 

0.0601 

0.0757 

0.0944 

0.0520 

0.0685 

0.0851 


T n 


0.0897 

0.3677 

0.7290 

0.0645 

0.1027 

0.1410 

0.0583 

0.0833 

0.1112 


r 


0.0826 

0.3618 

0.7118 

0.0585 

0.0961 

0.1333 

0.0541 

0.0766 

0.1037 


A F 


0.0781 

0.3258 

0.6809 

0.059 

0.0637 

0.0358 

0.0509 

0.0519 

0.0391 

7 

S(F, 

Ex) 

0.1117 

0.556 

0.9309 

0.0686 

0.1292 

0.1854 

0.0617 

0.0927 

0.1366 


Jo.9 


0.0784 

0.2153 

0.4461 

0.0549 

0.0812 

0.1032 

0.0529 

0.0666 

0.0881 


T n 


0.1084 

0.5255 

0.9181 

0.0664 

0.1266 

0.1846 

0.0608 

0.0906 

0.1314 


r 


0.1100 

0.5449 

0.9176 

0.0706 

0.1260 

0.1796 

0.063 

0.0974 

0.1395 


A F 


0.0988 

0.4600 

0.8758 

0.0602 

0.0806 

0.0668 

0.0533 

0.0612 

0.0556 

9 

S(F, 

Ex) 

0.1294 

0.7008 

0.9818 

0.0732 

0.1497 

0.2131 

0.0645 

0.1087 

0.1562 


Jo.9 


0.0857 

0.311 

0.6371 

0.0622 

0.0895 

0.1222 

0.0569 

0.0797 

0.1043 


T n 


0.1227 

0.6705 

0.9763 

0.0741 

0.1469 

0.2087 

0.0658 

0.1064 

0.1494 


r 


0.1284 

0.6798 

0.9708 

0.0697 

0.1364 

0.1914 

0.0679 

0.1072 

0.1520 


A F 


0.1068 

0.5828 

0.9476 

0.0630 

0.0952 

0.1032 

0.0557 

0.0710 

0.0804 

11 

S(F, 

Ex) 

0.1440 

0.8038 

0.9967 

0.0767 

0.1773 

0.2872 

0.066 

0.1285 

0.1994 


Jo.9 


0.0864 

0.3779 

0.7468 

0.0632 

0.1094 

0.1447 

0.0605 

0.0795 

0.1154 


T n 


0.1345 

0.7704 

0.9953 

0.0770 

0.1774 

0.2830 

0.0638 

0.1227 

0.1910 


r 


0.1419 

0.7775 

0.9931 

0.0717 

0.1563 

0.2494 

0.0657 

0.1202 

0.1844 


A F 


0.1204 

0.6867 

0.9857 

0.0695 

0.1193 

0.1452 

0.0571 

0.0849 

0.1033 

13 

S(F, 

Ex) 

0.1629 

0.8777 

0.9989 

0.0882 

0.2042 

0.3229 

0.0716 

0.1303 

0.2333 


Jo.9 


0.1011 

0.447 

0.8511 

0.0737 

0.1209 

0.1619 

0.0589 

0.0875 

0.1269 


T n 


0.1526 

0.8512 

0.9985 

0.0876 

0.2026 

0.3213 

0.0718 

0.1247 

0.2274 


r 


0.1557 

0.8494 

0.997 

0.0787 

0.1767 

0.2771 

0.0676 

0.1215 

0.2076 


A F 


0.1326 

0.7673 

0.9942 

0.0727 

0.1365 

0.1750 

0.0629 

0.0903 

0.1278 

15 

S{F, 

Ex) 

0.1764 

0.925 

0.9999 

0.0926 

0.2246 

0.3760 

0.0686 

0.1493 

0.2554 


Jo.9 


0.1074 

0.5271 

0.9143 

0.0705 

0.1259 

0.1802 

0.0613 

0.0895 

0.1430 


Tn 


0.1665 

0.9009 

0.9998 

0.0942 

0.2230 

0.3734 

0.0687 

0.1422 

0.2511 


r 


0.1800 

0.9067 

0.9995 

0.0862 

0.1903 

0.3217 

0.0664 

0.1362 

0.2345 


A F 


0.1439 

0.8261 

0.9984 

0.0767 

0.1466 

0.2085 

0.0606 

0.1005 

0.1474 
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4.2 The Two-Sample 


As mentioned in the introduction, Izadi and Khaledi ( 2012 ) proposed and studied a test 
for the two-sample problem based on kernel density estimation for testing Hq : F =* G 
against Hi : F <* G & F 7 C G. Their test statistic is 


A (F,G) = 


n 2 a, 


i= 1 3=1 


Xi - X , 


m 2 b r 


EE r <‘( 

i =1 3=1 


K 




where k is a known symmetric and bounded density function and a n and b m are two 
sequences of positive real numbers, k and a n are known as kernel and bandwidth, respec¬ 
tively. 

In this section, we compare the empirical power of our new test with the Izadi and 
Khaledi’s test when the kernel, k, is the density function of the standard normal distribu¬ 
tion and a n = n ~ 2 / 5 and b m = m ~ 2 / 5 . We know that the gamma and Weibull family are 
decreasing with respect to the shape parameter in the more IFRA order (cf. Marshal and 
Olkin, 2007, Chapter 9). Also, Izadi and Khaledi (2012) showed that the beta family with 
density function 

T a ~ 1(1 _ r \ b — 1 

Sip) = - - p( a ^ , a? e [0,1] , a,b> 0 . (4.15) 

is increasing with respect to b in the more IFRA order. So, to evaluate the power of 
the tests we use the gamma, Weibull and beta families denoted by G(a,/3), W{a,/3 ) and 
B(a, b), respectively, in the alternative hypothesis. In Table | 6 j we generated 10000 samples 
with sizes n = m = 20,30,40,50,100 from distribution F and G given in the table. We 
observe that the empirical power of our new test is greater than the empirical power of 
Izadi and Khaledi’s test when F and G belong to Weibull family and is smaller when F 
and G belong to the gamma and beta families. So, our new test is comparable to Izadi 
and Khaledi’s test. 


5 An application 

In this section we apply our test on a data set from Nelson (1982, page 529) which is 
a life test to compare two different (old and new) snubber designs. Let F (G) be the 
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Table 6: The empirical power of the tests 5(F, G ) and A (F, G ) 



Distribution 




n = m 



F 

G 

test 

20 

30 

40 

50 

100 

G{ 3,1) 

G(1.5,1) 

5(5, g) 

0.415 

0.582 

0.6764 

0.7692 

0.9589 



A (5,G) 

0.470 

0.642 

0.7554 

0.827 

0.9712 

G(4,l) 

G( 2,1) 

5(5, G) 

0.435 

0.570 

0.681 

0.785 

0.966 



A (5,G) 

0.492 

0.672 

0.755 

0.825 

0.968 

W (3,1) 

W(1.5,1) 

5(5, G) 

0.7946 

0.9346 

0.9804 

0.9940 

1 



A(5,G) 

0.6714 

0.8562 

0.9446 

0.9796 

1 

W( 4,1) 

W(2,l) 

5(5, G) 

0.811 

0.9312 

0.9766 

0.993 

1 



A (5, G) 

0.72 

0.893 

0.954 

0.985 

1 

B{ 1,1.5) 

5(1,3) 

5(5, G) 

0.1292 

0.1736 

0.2332 

0.2692 

0.4396 



A (5, G) 

0.1452 

0.2048 

0.2806 

0.3314 

0.5504 

5(1.5, 2) 

5(1.5, 5) 

5(5, G) 

0.166 

0.209 

0.278 

0.364 

0.597 



A (5, G) 

0.383 

0.517 

0.585 

0.653 

0.938 
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distribution of lifetime old (new) design population. In Fig. 5, Izadi and Khaledi (2012) 
plotted TTT-plots for both data sets of old and new design. The graphs anticipated IFRA 
populations for both populations. 

Now we apply our IFRA test on the two data sets. Using our one sample test, we 
get that \/12 n{8p — 3/2) = —3.579222 and y/12m(Sc — 3/2) = —3.085525 which are less 
than —2.376441 (the critical value at level of significance a = 0.01 from Table [[]). So, 
our test reject exponentiality of both population in favor of IFRA-ness. To compare two 
populations with respect to more IFRA order, the test statistic value of the two sample 
problem is \/N5(F,G) / a F ,G = —0.3762 ^ —2.326348 = £ 0 .oi- So, at level of significance 
a = 0.01, the equality of two populations in more IFRA order is not rejected. 

6 Summary and Conclusion 

In order to evaluate the performance of an appliance, we need to compare its aging be¬ 
havior with some distributions such as exponential, Weibull, gamma, linear failure rate 
distributions. The notion of the star order ( denoted by <*) is one of the useful tools for 
this comparison between two distributions. 

In this paper, we have introduced a new simple test for the problem of testing H 0 : 
F =* G against Hi : F <* G and F 7 U G. 

In the one-sample problem, let Xi ,..., X n be a random sample from F and G = Go 
where Go is a known distribution. Hq is rejected at level of significance a, for large sample 
size, if — bG 0 )/&G 0 < z ai where 

.. EE max {Ah, Xjj 

r _ E g 0 [max{Ad,X 2 }] i _ 

° G ° ~ -’ ° F ~ - 7 - Hv - 

f’Go n[n - 1)A 

and 

i f o S 

<4 = — x Var Go XG 0 (X) + / tdG 0 (t) - ^X 
^G 0 \ J X z 

In particular, when Go is an exponential distribution, the null hypothesis in favor of IFRA- 

ness is rejected, if \/l2n(S F — 3/2) < z a . The exact null distribution of the test statistic 

has been obtained and, for small sample sizes 2(1)40, the exact critical points of the 



16 








test statistics have been computed. Based on Pitman’s asymptotic relative efficiency and 
simulated power, we have compared our test with the tests given by Deshpande (1983), 
Kochar (1985), Link (1989) and Ahmad (2000). The results showed that our test relatively 
dominates the other tests. 

In two-sample problem, let AR,..., X n and l'i,..., Y m be two random samples from F 
and G respectively. For large sample sizes, we reject H 0 in favor of H i if 

VN{5f — &g) /&F,G < A* 

where N = n + m and &f,g has been given in (13.101) . Using simulation study, we have 
shown that our test in this case is comparable with the test of Izadi and Khaledi (2012). 
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